INTERPOLANTS OF LATTICE FUNCTIONS FOR THE 
ANALYSIS OF ATOMISTIC/CONTINUUM MULTISCALE 

METHODS 



C. ORTNER AND A. V. SHAPEEV 

Abstract. We introduce a general class of (quasi-)interpolants of functions de- 
fined on a Bravais lattice, and establish several technical results for these in- 
terpolants that are crucial ingredients in the analysis of atomistic models and 
atomistic/continuum multi-scale methods. 



1. Introduction 

In this report, we describe several classes of continuous interpolants or quasi- 
interpolants of lattice functions that have proven useful in recent and ongoing work 
[HI HOI [11] for the analysis of atomistic models of crystalline solids and for for the 
construction and analysis of atomistic/continuum multiscale methods [21 [13], [3 [6j [12]. 

A key ingredient that made the analyses in [8] [TOj CHI El E] possible is the con- 
struction of continuous representations of discrete object. A novel recent idea that 
allowed substantial progress in the direction of such constructions is the bond density 
lemma [12] , which allows explicitly compute the "density of bonds" in a certain sub- 
sets of M. d . The lattice interpolants we define and analyze in this paper are designed 
to go hand in hand with these constructions. 

Although many results related to the ones we present here can be found in the 
numerical analysis literature, we have not found the precise statements we require 
or the generality that we seek (as a special case our analysis covers multi-linear and 
multi-cubic splines for which many of our results are standard). Hence, we present 
a complete derivation of all results. Moreover, our function space setting seems new 
as well. 

2. Interpolation of Lattice Functions 

We fix a space dimension d G N and a range dimension m G N. The goal of 
this section, and core of this paper, is to introduce spaces of discrete and continous 
functions, respectively, on the continous domain M. d and the discrete domain 7* d , and 
to provide tools to transition between these two classes. 

2.1. Interpolants of lattice functions. We denote the set of all vector- valued 
lattice functions by 

9/ := {v : Z d -»■ W n }. 
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To facilitate the transition between continuous and discrete maps we introduce two 
(quasi-)interpolants of lattice functions, and will later introduce a third nodal inter- 
polant. 

Our starting point is a nodal basis function ( G W /1 '°°(]R d ; M) associated with the 
origin; that is ((0) = 1 and = for all £ G Z rf \ {0}. There is considerable 
freedom in the choice of (, however, certain choices are particularly natural; see 
Section |2.2| The nodal basis function associated with £ ^ is the shifted function 



£(• — £), which yields the interpolant 

v(x):=^2v(£)C(x-(i), forveW. (1) 



Next, we define a quasi-interpolant obtained through convolution of v with 

v(x) :— (C * v)(x) — / ^(x — x')v(x') dx' . (2) 



We call v a quasi-interpolant since, in general, v (£) ^ v (£). This can be seen from 
the alternative definition 



v\x) 



5>(0C(:c-0> where CW := (C * W - (3) 



For future reference we define the support sets 

w$ :=supp(C(»-f))> and ^ := su Pp(C( # -£))■ 

Throughout, we make the following standing assumptions on £: 

(Zl) Regularity: ( G I^ 1 ' 00 

(Z2) Locality: (and hence Of) is compact 

(Z3) Order: C{ x — ( a + ^ ' £) = a + • £ for all a G IR, 6 G R d . 

(Z4) Invertibility: ((0) = 1 and C(0 = for all ( G Z d \ {0}. 

Remark 1. 1. The assumption (Z3) ensures that interpolants of affine functions 
are again affine: if = a + b-£, then -u(x) = a + b-x. When applied to continuous 
functions in § |4| this ensures that this interpolant is second-order accurate. There, 
we will also see that a nodal interpolant from the space {v} is third-order accurate. 

2. The assumption (Z4) is necessary not only to ensure that v interpolates v, 
but also to ensure that the operation v 1— > v is invertible. Indeed, the extended hat 
function in ID 

1/3 -1<x<1 
l/3(2 + :c) -2<x<l 
1/3(2 -x) l<x<2 
Id > 2 



satisfies (Z1)-(Z3), yet for the function u = (— 1,0, l) per (periodic repetition of 
(-1,0, 1)) we have u = 0. □ 

Lemma 1. Let v G <% , then v G W^(R d ; R m ) and v G W^(R d ; R m ). 

Proof. The basis function ( has compact support and belongs to W 1,oa . Hence, in 
every compact subset of lR d , v is a finite linear combination of Lipschitz functions, 
which implies that v G . 
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Figure 1 . Subdivision of the unit cube leading to a regular simpli- 
cial partition of M. d that is translation invariant and symmetric about 
the origin. 



A straightforward computation shows that 

V 2 C(x) = J VC(x - x) ® VC(x') dx, 

hence V 2 C is Lipschitz continuous and in particular ( G W 3,OQ with compact support. 
The same argument as above applies to prove that v G W£ c °°. □ 

2.2. Examples. 

(1) Ql-interpolation: possibly the most natural choice for the basis function ( 
is the tensor product 

d 

C(x) := Jjmax(0,l-|zi|). (4) 

i=i 

In this case, supp(C) = [—1, l] d is compact and ( is piecewise multi-linear; 
{v | v G ^ } is the space of tensor product linear B-splines on the grid 
and {v | v G % } is the space of cubic tensor product B-splines on the grid 
(see |U). 

(2) Pl-interpolant: In some cases, such as in [llj, it is convenient if Vu is 
piecewise constant. This is also possible within our framework. Let ST be a 
regular simplicial partition of IR d , which is translation invariant (£ + = ^ 
for all C, G 7j d ) and symmetric about the origin {—2T = 3?), and let ( denote 
the Pl-nodal basis function with respect to this partition. Then it is easy to 
see that ( satisfies (Z1)-(Z4). 

For d — 1,2 such sub-divisions are straightforward to construct. For d = 3, 
the subdivision of the cube shown in Figure [T] yields a partition with the 
desired properties. 

2.3. Basic stability results. In the following theorem we collect all basic stability 
results for the first two lattice interpolants. In these statements the norms may take 
the value +oo. 

Theorem 2. There exist constants Co,Ci,C2 > ; independent of p, such that 
CqCi\\u\\p < Ci\\u\\ep < \\u\\lp < \\u\\lp < c 2||w||^p Vw G pG[l,oo). (5) 
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Moreover, there exists a constant c\ > 0, independent of p, such that 

c'i\\Vu\\ L p {Rd) < ||V£t|| LP(Kd) < ||V«|| L p (R «i) VmG^, pe [l,oo). (6) 
All inequalities in ^ and (J7|, except the first in each display, also hold for p = 00. 
We establish this result in a series of lemmas. 

Lemma 3. Let ^ be a finite- dimensional vector space and let Si,S 2 '■ ^ — > K 
be semi-norms. If ker(S 2 ) C ker(Si) then there exists a constant C > such that 
Si(u) < CS 2 {u) for alined. 

Lemma 4. There exist constants Ci,c 2 > such that 

Cl||w||#> < ||u||z,p < ||u||z,p < C2 1 1 1i 1 1 £p 

Proof. 1. Proof of \\u\[u> < \\u\\lp: This statement follows simply from the fact 
that u = ( * u, where ( > and J (dx = 1 (see, e.g., (3], Sec. C.4]). 

2. Proof of \\u\\ LP < c 2 \\u\\ P . Fix £ E Z d and let u := £ + [-1/2, l/2} d . Let 
Si(u) := \u(£)\ and S 2 (u) := \\u\\jj>( u -), then Lemma [3] implies that 

< c 2 ||u||^. 

An application of Holder's inequality yields 




Since ( has compact support, we can bound 

(Ek>-oi) p/p '<^, 

undidendently of x and p. Hence, we obtain 

\\u\\% <C [ T\<0\ P \((x-0\dx = C\\u\\ £ 4(\\ L i. 

JRd fez* 

3. Proof ofc l \\u\\iP < \\u\\lp: If m = in £ + [-1/2, l/2] d then u(£) = 0, hence 
Lemma [3] implies that 

\u(£,)\ < C , ||'u|| i i( ?+ [_ 1 /2,l/2] d ) < C||w||LP(g+[-l/2,l/2] d )? 

where C is independent of p. Summing over £ e Z d (and with suitable modification 
for p = 00) yields 

\\u\\iP < C\\u\\lp, for u G □ 

It remains to prove that co||m||#> < ||u||^p. We will first prove this for p = 2 
and then use Fourier analysis to extend the result to general p. It is convenient 
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to interpret this analysis in terms of the operator ^ : £ p (Z d ;W m ) — > £ p (Z d ;M. m ), 

Lemma 5. ^ : £ 2 — )■ £ 2 is an isomorphism; that is, there exists Cq > such that 

c olMI^ 2 < ll^ll^ 2 ^ u ^ ^ ■ 

Proof. Let u G £ 2 . Using the same argument as in step 3 of the proof of Lemma E] 
we obtain 

\\u\\p < c^WuWv, 

hence, we can estimate 




=c o 1 ^2 u (o -«(o ^ c o 1 ii^iU 2 ii^iU 2 - □ 

Lemma 6. For any p G [l,oo) ; ^ : £ p — > £ p is an isomorphism; in particular, 
there exists cq > such that 

Co\\u\\ep < 11^11^ Vn G 
Proof. We use the representation 

where m(£ — rf) := f — — x) dx. 

Next, we introduce a discrete Fourier transform : £ 2 — > L 2 ([0, l] d ), 

With this definition = ,^~ x {m^u\ where rh = J^~m G L 2 ([0, is the Fourier 
multiplier associated with ^ . Since is an isomorphism between £ 2 and L 2 ([0, l] d ), 
we can express 

W^Wp = (minm(a))- 1 , 

a 

and from Lemma[5]we can deduce that min Q m(a) > c > 0. Since m is analytic and 
bounded below it follows that g(a) := (m(a)) _1 is an analytic function on [0, l] d . 
Next, we show that g G i 1 , where g := ,^~ x g. Let k > d/2, then 

ii* < Il(i + \er k L ll(! + KHVOIU 

< H(l + l£lT1^ ||(id-A)^|| LW) . 

The first term is bounded since k > d/2, while the second term is bounded due to 
analyticity of g. 

Since g G £ x it is easy to see that the operator £f : £ p — >■ £ p , 

(Sf*)(0:=5>(£ -*?)«(»?)> 
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is well-defined, for all p £ [1, oo]: For p £ {1, 00} it is well-known that \\&\\&> = \\g\\p 
(where \\&\\ep is the £ p -operator norm), and hence, by the Riesz-Thorin interpolation 
theorem, ||£f||#> < \\g\\p for all p £ [l,oo]. 

It only remains to show that £f = c io~ 1 . This follows from the fact that, by 
construction, £f = on £ 2 and density of functions with compact support in t v 
for p < 00. □ 

Lemmas |3]-[6] establishes ([5]). We will now prove the corresponding result for the 
gradients. 

Lemma 7. There exists a constant c[ > 0, independent of p, such that 

ci||VM|| z , < ||Vu||lp < \\Vu\\ LP VueW, pe[l,oo). (7) 
The second bound also holds forp = 00. 

Proof. 1. Upper bound. Since u = ( * u it follows that Vm = C, * Vu. Hence, the 
upper bound follows, e.g., from P, Sec. C.4]. 

2. Lower bound. Let v £ % be defined by v (£) :— «(£), then 

IIVuHlp < ||Vw - Vv\\ LP + ||Vv||lp. (8) 

2.1 Estimating ||Vw||lp. This term can be estimated by a simple semi-norm 
equivalence argument. Fix a hypercube Q £ J2 and define Si(u) := ||Vw||£j>(q) 
and 5*2(u) := ||Vw||lp(q). Si, S2 are seminorms on % and involve only a finite 
number of degrees of freedom. Moreover, if Si{u) = 0, then u is constant in Q, and 
consequently, Si(u) = 0. According to Lemma |3j there exists a constant C > 0, 
such that Si(u) < C 82(11) for all u £ y . Due to translation invariance, the constant 
must be independent of the choice of Q. Summing over all cubes we obtain 

||Vv|| L p (R d) < C||Vu|| LP ( K d). (9) 

2.2 Estimating ||Vu — Vw||lp. Suppose, first, that 1 < p < 00, then 

HV«Cp= / |5>(£)VC0r-£) P clx 

jRd ' fGZ- 

< / (E hoi p ivc> - 01) ( E i v ^> - ^)i) P/P ' dx - 

Since |VC| < 1, and since the number of lattice sites £ such that |VC(x — £)| > is 
bounded independently of s £ IR d , we obtain 



Vi0||k < Ci E H0\ P I |VC> - 01 da;. 



Using again < 1 and the compact support of VC we obtain 

II Vw||ip( K d) < C 2 ||w||^(Zd), 

for some constant C2 that is independent of w. Applying ^ and reinserting the 
definition of w, we arrive at 

||Vw - Vv|| iP ( R d) < C 3 \\u - v\\ip(z,d). (10) 
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Since v is a quasi-interpolant of u we can bound the right-hand side of (10) in 
terms of Vu. For a fixed lattice site £ G 7j d , we now obtain 



m-m = m-m= / [m-^)]c(x-odx. 

Since f Rd ((x — £) dx = 1, this implies 

|fi(0-«(0| < ( / \m-Hx)\ P C(x-0dx) 1/P =:S 1 (u). (11) 

V ./rod / 



Moreover, let 



( [ \Vu\ p ((x-0dx 



VP 

) 



then Si, S2 are both semi- norms on J^. Clearly, if S2(w) = then u is constant in 
which implies that also Si(u) = 0. Applying Lemma [3J we obtain a constant 
C 4 > such that 

Si(u) < C 4 S 2 (u) Vm G 

Due to translation invariance, the constant C4 is independent of £. Combining this 
result with (11), (10) and the fact that £(• — £), £ G Z d is a partition of unity, we 
obtain 



\vu-vv\\ p LP <c p 3 j2Ho-m\ 



< Cf CJ £ / \Vu\ p C(x -t)dx= (C 3 C A y\\Vu\\ p LP d (12) 

This concludes the estimate of || Vw — Vu||jy». 

Combining (J8J), (pi), and (12) yields the lower bound in Q, for p G (l,oo). For 
p — 1 a minor variation of the argument gives the same result. □ 

Finally, we state some useful embedding results. 

Lemma 8 (Embeddings). (i) Let p, q G [1, 00] and i < j G {0, . . . , 3}; then there 
exists a constant C such that 

\\V j u\\ L p {Q) < CWV'uWl^q) V«6^, VQeQQ. (13) 

(ii) Moreover, if p> q, then 

\\V j u\\ L pmd) < CUV^Hi^Rd) Vu G W '. (14) 



Proof, (i) The first estimate results from local norm-equivalence in each element 
Q G £}, due to the fact that the spaces {u\q \ u G % } are finite-dimensional and 
Q- independent. 

(ii) The second estimate is a consequence of the embedding l q C l v . □ 

2.4. Discrete Sobolev spaces. We define the discrete Sobolev spaces W 1,p ', p G 
[l,oo], by 

:={«g^|mG W l ' p ), 

equipped with the norm ||w||^i, P := In view of the embedding (14) the 

discrete Sobolev norm is equivalent to the £ p -norm, and hence, W 1,p = l p . We will 
introduce the more natural homogeneous Sobolev spaces in the next sub-section. 
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2.5. Discrete homogeneous Sobolev spaces. We introduce function spaces that 
naturally arise in the analysis of atomistic models and their approximations. We 
define the (semi-)norm 

||m||^i,p := ||Vm||lp, for u G W , pe[l,oo]. 

Since || • \\wi,p does not penalize translations, we define the equivalence classes 

[u] ■= {u + t \t G R d }, ioiueW. 

We will not make the distinction between u and [u], whenever it is possible to do 
so without confusion, for example, when all quantities involved in a statement are 
translation invariant. 

For p G [l,oo] we define the discrete function space 

^ x >f : = {[u] \ u G \\Vu\\ L p < +00}, (15) 

which we analyze in the following results. We will also require the space of displace- 
ments for which the displacement gradient has compact support, 

^0 := {H £ ^ supp(V-u) is compact}. 

Proposition 9. ^ 1,p is a Banach space. For p G [l,oo), is dense in ^ 1,p . 

Proof. The proof for continuous homogeneous Sobolev spaces [9J Prop. 2.1 and 
Thm. 2.1] translates verbatim. 

In these results, density of (equivalence classes of) displacements with compact 
support was shown, which is a small class than & . Hence the case d = p = 1 
was excluded. Here, we admit the slightly larger class and it is straightforward to 
show that p = d = 1 now also included in the result. (For d > 1 the two classes 
coincide.) □ 

Next, we analyze the convolution operator ^ on % l ' p . 

Lemma 10. Let p G [l,oo) then & : Of 1 * W l > p , tfu : = (u(0)?ez d ; iS an 
isomorphism. Moreover, we have the stronger estimate 

ll^u-nll <C||Vn|| LP . (16) 



Proof. Fix u G < W l ' p . We need to show that there exists u* G % 1,p such that u* = u. 
This is equivalent to the statement that there exists w G such that w = u — u. 
By repeating the argument following ( 11 ) it follows that, in fact, u — u G £ p . Since we 
know that ^ : i p — > i p is an isomorphism, we can define w := "if — u). Applying 

it follows that <*f : W 1 * -> ty 1 * 

□ 



Lemma^ we obtain (16). Since ||Vw||lp < C\\w\ 



is indeed an isomorphism. 



2.6. A smooth nodal interpolant. In view of Lemma [6] and Lemma 10 
now define the nodal interpolant 

Iu:=iPy forwG^UF, pe[l,oo). 



we can 



(17) 



Lemma 10 immediately implies the following norm-equivalence result. 



Lemma 11. There exist constants Cq,ci, such that 

c ||Vm|| L p < \\VIu\\ LP < cAVuhp Vm G ^ 1,p , 



p e [1,00). 
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3. Discrete Deformation and Displacement Spaces 

The goal of this short section is to make rigorous the meaning of the far-field 
boundary condition for a discrete deformation field y : Z d — > M d , 

y(0~Af as|£|->oo. (19) 



In order to rigorously define the far-field boundary condition (19) it is useful to 
analyze the asymptotic bahaviour of displacements u G ^ 1,p . 

Proposition 12. Let p G [1, oo) and u G <%f l,p . 

(i) Ifp > d, then \u{£) -«(0)| < C||«||*i*|^| 1 /p'; 

(ii) Ifp = d, then \u(l) - u(0)| < cjjujj^^ + log |£|); 

(iii) If p < d, then there exists a uq G [u] such that uq G £ p * . 

Proof. First, we note that < W 1)P C ^ 1,0 ° for all p G [l,oo). The result therefore 
follows directly from |9, Thm. 2.2]. □ 

Proposition [9] allows us to give a clear interpretation to the far-field boundary 
condition (ph. Let A G M+ Xd , y A (x) := Ax for x G M d , and let 



— > 0, or equivalently, ?/(£) = A£ + o(|£|), as |£| — > oo. 



^^{MlyG^Jy-^G^ 1 ^}. (20) 
From Proposition |9] we infer that, for all y G [y] G ^ 1,p , 

|y(0-Ae| 

lei 

Motivated by this discussion we may confidently take as suitable classes of 
admissible deformations, or equivalently, fy? 1,p as suitable classes of admissible dis- 
placements. 

4. Approximation Error Estimates 

4.1. Nodal interpolants. In this section we prove two useful interpolation error 
estimates. We define the nodal first-order interpolant and extend the nodal third- 
order interpolant, for v G C(IR d ; IR m ), by 

Iv[x) := ^ v (0(( x ~ 0> and M x ) := (H v \z d ))( x )i 

where the latter is well-defined provided that v\%d G % 1 ' p ) p G [1, 00]. Moreover, for 
v = yA + u, u G we can define Iv = y& + Iu. 

Lemma 13. Let k G {1, 2}, and p G (d/k, 00] if d > 1; then there exists a constant 
Cj such that, for all v G W k ' p (Q), and Q G B, 

\\Iv - v\\ W k, P{Q) < C]\\V k v\\ L p {Q) . 

In particular, if v G W^f with V k v G L p , then 

\\Iv — u||jyfe,j> < Cj\\ \7 k v\\ LP . 
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(ii) Let k G {1,2,3,4}, and p G (d/k, oo] if d > 1; t/ien t/iere exists a constant Cf 
such that, for all v G W^' c p witt G W^" 1 '*', [w| 2 d] G ^ 1>p and 



\Iv — u||v^fe,P < Cj\\ V k v\\ LP . 



Proof, (i) The first statement follows from a standard Bramble-Hilbert argument; 
the condition p > d/k (p > 1 if d — 1) ensures that the nodal values are well-defined 
and hence the interpolation operator is stable; see [1] for various examples of the 
argument. 

(11.1) We first prove that [v\%d] G e % 1,p . From (i) we know that Iv — v G W k,p and 
in particular V/f G L p . By definition, this implies that [v \%d] G ^ 1,p . Hence, Iv is 
well-defined. 

(11.2) The interpolation error estimate for / requires some care since the inter- 
polation operator I and hence I are defined through a linear system, i.e., they are 
non-local. To prove this result we choose an arbirary w G ( W X,P , and estimate 

\\Iv — V \\w k 'P — \\I V — + \\w — v \\w k .P- 



Clearly, Iw = w. Using (14) we obtain 

\\Iv -w\\ W k, P = \\I(v - w)\\ W k, P < C\\I(v - w)\\ LP < C\\I(v - w)\\ LP , 

where the last inequality follows from Lemma [6] and from Lemma [2] The assumption 
that p > d/k if d > 1 ensures that / is stable, that is, \\I(v — w)\\lp < C||t> — w\\ w k <p- 
Thus, we have obtained that 

\\Iv — v\\ W k, P < C\\w — V \\Y/k, P . 

Now choosing w to be a suitable quasi-interpolant gives the desired result; this is 



provided by Lemma 17 in the next sub-section. □ 



Remark 2. Let p > d/k if d > 1 and arbitrary otherwise (e.g., p = 2 and 
d G {1,2,3}). Let u G e % l > p be a discrete displacement, e.g., the equilibrium dis- 



placement of an atomistic model. Then Lemma 13 (ii) implies that, for any function 
u G W k,p interpolating u in lattice sites, we have 

\\V k Iu\\ L p < C||V fc n|| LP , 

for some generic constant C . Thus, up to a generic constant, Iu is the "smoothest" 
interpolant of u. In particular, V k Iu is a canonical measure for the smoothness of 
the discrete map u. □ 

4.2. A quasi-interpolant. We are left to define and analyze the quasi-interpolant 



used in the proof of Lemma 13 (ii). The resulting interpolant will also be interesting 
in its own right since it will remove the restriction p > d/k imposed on the nodal 
interpolant /. 

The main ideas of our construction are standard, however, the details require 
some care. The idea, following Clement [2], is to construct a bi-orthogonal basis 
function with compact support such that 

?MCtt-*)d* = { J; « t =° wise , (2i) 
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and to define 

J«(a:):=X)(C**«)(OC(x-0- ( 22 ) 



Lemma 14. There exists (* G L°°(R d ) with support supp (* = Uq, satisfying (21). 

Proof. Let X := {£ G Z d | |u>^ fl u \ > 0}, and let S£ := span{^ | £ G X}, where 
= CO — 0- We seek (* of the form 



C(x) := 



0, otherwise, 



for some parameters a^. 

Since the functions Q are linearly independent in M. d , ( cannot be written as a 
linear combination of G X. Therefore, the mapping 

Sex 

is well-defined and a linear functional on 3£ . Thus, by the Riesz representation 
theorem, for equipped with the L? inner product, there exist coefficients (a^)^ e x 
such that 

'1, i = 0, 



0, otherwise. 



Since supp(C*) = Wo the result follows. □ 

Lemma 15. Let Ci arid C2 fre two different functions satisfying (Z1)-(Z3), let 
Ci := Ci * Ci> an d kt p(x) be a cubic polynomial. Then 

£K0 -&(*-£)). 



Proof. Denote 5( := (2 — Ci- From (Z3) we have that 

E set* -0 = 0. 

Thus, there exist Lipschitz functions /& : lR d — >■ lR d , = 1, . . . , d, with compact 
support, such that 

d 

k=l 

where Dku(x) := u{x + e&) — 

Further applying summation by parts, we have that, for all £ = l,...,d and 



x G R d , 



d 



0= E^^-^E^E^^-o 



it=i 



E - = EEw w (^ - = E - o- 

fc=i £ £ z d fe=i 
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From the last identity, we have that 

f^\x) = j^D k f^\x) 
k=i 

and hence 



5((x) = E P,Ihf IJ - i 



i,k) t 
k,e=i 

We can now return to proving the statement of the lemma. Using the represen- 
tation 

Ca = Ci * Ci + 2Ci * S( + $( * S( = Ci + (2Ci + *C) * K 



we have 



Eko(c 2 (^-o-0-o) 



= E E P(0 (( 2 Ci + *C) * D t D k fM) (x - o 

= E E p(0^*((2Ci + ^) * f m )(x - o 

£ G zd e,k=i 

d 

= E E (^*p(0) ((26 + 50 * / (Afe) ) (x - o- 

§ez d £,fc=i 

Since D e D k p(^) is a linear function for all £ and fc and Ci satisfies (Z3), we have that 

E E (^*p(0) (2G * / (£>Jfc) ) (x - = 2 E / D i D #(* - y)f m (y) ^ 

£ & zd e,k=i e,k=i J 

is indeed a linear function of x. It can be shown that the remaining part is zero: 

E E(^^(0)K*/ ( '' fe) )(^-o 



= EE t^^^toKP^/^j^^o 

^ e ^ d ^fc=i^,fe'=i 

since D e D k >D e D k p(£) = 0. □ 

Lemma 16. Lei fre a cnfrzc polynomial, then there exists w G ^ snc/i £/ia£ 
w — p. 
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Proof. Assume, without loss of generality, that m = 1. 

First, suppose, that p(x) = a T x + c is linear. Let v(£) := a T £ + c, then by (Z3), 
v(x) = p. Hence ti = (*?; = pas well, since the convolution with ( preserves linear 
functions. 

Now notice that, if £ is the Ql-nodal basis function given by (|4]), then {w \ w G % }, 
is the set of all multi-cubic B-splines and therefore contains p [3]; i.e., p = v for some 
wGf. For a general (, Lemma 15 guarantees that v — p is a linear function, where 
v G % is the function we just constructed. Hence letting w := v — (v — p)|z<» and 
using again the fact that convolution with ( leaves linear functions invariant, yields 
w — v — (v — p) = p. □ 



Theorem 17. Let J be defined by (22) with £* from Lemma 14 and let < j < 



k < 4; then, for all v G Wjf , G W*- 1 *, we /jaw 



I LP ^ ^ 



|V fc 



D if. 



Proof. For any function u; = X^ez d U '(C)C( ;E — 0) the definition of J guarantees that 
Jw = w, that is, J is a projector. 

We begin by applying a triangle inequality. For any Q 6 J, < j < k < 4 and 
w, w G , we have 

||V J '(J« - < l|V J (J^ - w)\\lp{Q) + II V> - v)\\ LP{Q) . (23) 



Using the inverse estimate (13), and a straightforward computation (or, a local 
norm-equivalence argument) we obtain 

\\V J {Jv-w)\\ LP{Q) = \\V J J{v-w)\\ LP{Q) < C\\J{v-w)\\ LP{Q) < C\\J(v-w) WtpptraQ, 
where oj'q := U§ez d nQ ^5- For eacri £ G Z d fl cDq, we have 



\J(v-w)(0\ 



C( x — 0( v ~ w)(x) dx 

oj 5 



<C\\v-w\\ LP{Cl6) . (24) 



Combining (|24h with (23), we obtain 



||V J (Jf - U)||£P(Q) < C(\\v - w\\ L v{u Q ) + \\^(v - w)\\ LPiQ) ), 

where Cjq D Q'q is a compact set of uniformly bounded diameter. According to 
Lemma 16 we may choose w to be any cubic polynomial. Minimizing the right-hand 
side over all cubics yields 

\\V^Jv-v)\\ LP(Q) <C\\V k v\\ LP{C0Q) . 

Summing over all Q G J2 and estimating the overlap of the sets Oq completes the 
proof. □ 
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